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ABSTRACT 

Effective mechanical properties for large repetitive frame- 
like structures are derived using straight forward combinations 
of strength of material and orthogonal transformation techniques. 
Once the actual structure is identified symmetry considerations 
are used in order to identify its independent property constants. 
The actual values of these constants are constructed according to 
a building block format which is carried out in the three consec- 
utive steps: (a) All basic planar lattices are identifed (b) 

effective continuum properties are derived for each of these 
planar basic grids using matrix structural analysis methods and 
(c) orthogonal transformations are finally used to determine the 
contribution of each basic set to the overall effective continuum 
properties of the structure. 
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I. Introduction 


In racent two papers (1/21 we introduced a straightforward 
construction procedure in order to derive continuum equivalence 
of discrete pin jointed repetitive structures. Broadly speaking 
we outlined the method as follows: Once the actual structure 

was specified symmetry considerations were used in order to 
identify its independent property constants. The actual values 
of these constants were constructed in accordance with a build- 
ing block approach consisting of '-he following three consecutive 
steps: (a) all sets of parallel members were identified, 

(b) unidirectional "effective continuum" properties were de- 
rived for each of these sets and (c) orthogonal transformations 
were finally used to determine the contribution of each set to 
the overall effective continuum properties of the structure. 

Here the term properties is general and includes mechanical 
. tif fnesses) , thermal (coefficients of thermal expansions) and 
material densities. The method was then applied to a variety 
of structures. 

In the present paper we extend the analysis of [1,2] in order 
to derive the effective properties of rigid- jointed (fraue-like) 
repetitive structures. This differs substantially from the 
truss-like structures in that we here include the influence of 
inplane bending rigidities to the structure. The construction pro- 
cedure will differ in that the rod's unidirectional properties 
will no longer be adequate to derive the overall properties. 

The fact that the individual rod in a rigid- jointed array can 
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resist in plane bending dictates that the smallest suJt/-cell of 
the structure which vfill be used for the building block 
approach will no longer be unidirectional and thus have to be 
two-dimensional substructures. Here the most identifiable basic 
two dimensional frame structures are the (0°, 90®) and 
(0®, ^ 60®) lay ups. Effective properties for the sub-cells 
will be constructed using the direct analysis method which is 
also known by the matrix structural analysis method (see, for 
exairple (3-5]). This method, which uses simple and straight- 
forward strength of material techniques, constituties two- 
dimensional generalization of the one-dimensional area weighted 
properties approach of [1,2], The derived effective properties 
for such substructures will then be used in a bt Iding block 
format in order to derive the effective properties of more 
complicated two and three dimensional structures. This last 
step will be done by employing the orthogonal transformation. 

In summary we thus outline the procedure of constructing 
effective properties for frame-like repetitive structures as 
follows. Once the actual structure is identified symmetry 
considerations are used in order to identify its independent 
property constants. The actual values of these constants are 
constructed according to a building block format which is 
carried out in the three consecutive steps: (a) All basic 

planar lattices are identified (b) effective continuum 
properties are derived for each of these planar basic grids. 

Here a representative repeating cell is isolated and studied 



by the direct method noting that the effect of the joints’ 
rigidity is taken into consideration and (c) orthogonal trans- 
formations are finally used to determine the contribution of 
each basic set to the overall effective continuum properties 
of the structure. 

Since the inclusion of bending rigidities do not influence 
the thermal expansion of the structure, the thermal expansion 
properties derived in il,2] for the truss are identical to those 
of corresponding frame. Accordingly in what follows we con- 
centrate on deriving the elastic properties of tiie frame structure. 
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II. Orthoaonal Transformations 


As was pointed out earlier the actual values of the total 
structure's effective continuuir properties are determined from 
the individual contribution of each two-dimensional subset. The 
individual subsets contribution is obtained by u three-dimensional 
coordinate transformation. Before we proceed to describe the 
transformation, however, we shall first state the relevant stress- 
strain relations of elastic bodies. 

The stress-strain relations for a general linear elastic 
body are written in the compact ft rm 


'■"’ij ” ^ijke*'ki^ 


i,j,k,t = 1,2,3 


( 1 ) 


where c. . and c, . are 
1] k? 

tensors, respectively 
solid . 

^or future format 
in its expanded form 


tne components of the stress and strain 
and C. , are the stiffness tensor of the 

1 j k t 

references we shall rewrite equation (1) 
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Since C. , is a fourth-order tensor it obevs the transformation 

1 j K > 


[1,6,7 


ijkiL pqrs pi qj rk si 


where 


are components of the orthogonal transformation tensor which 
transforms the unprimed to the primed coordinates. Accordingly, 

I 

e. . is the cosine of the angle between the x. and the x. axis. 

The relation (3) hold equally well for either continuous or 
discrete structures. The numerical values of the aprropriate 
C. , entries will depend, however, upon the specific structure 
under consideration. Since we are interested in analyzing frame- 
type structures that are constructed from smaller subsets, it is 
expected that each s\ibset will contribute to its overall properties 
If a structure has n different subsets then equation (3) 


= 'Sqrs Saj > <5: 

Once the direction cosines of each subset are identified 
the sum over all of these subsets yield the final properties 


C. ; (C. ,, J 

ijkx. ijkjc 

m=l m 
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Ill . 


Basic Planar Grids 


We shall use the 'direct method" to find the properties of 
the equivalent continuum of two basic planar qrids. This 
approach is the reverse of that used by McCormick 18], 

McHenry, [9] and Hrennikoff [10], who describe a procedure 

for modeling problems in plane stress analysis with one dimensional 

elements . 

The main idea behind the direct method is to equate the 
displacements of the nodes of the model to the displacements of 
the corners of the continuum plate element under the same loading 
conditions. The sign convention for the displacement and stress 
resultants used in the present study are shown in sketch la,b. 
a) (0^, 90°) layup 

We consider a plane network which is formed from a large 

number of orthoaonai ly intersectino beams ricidly jointed at 

their intersections ae shown in figure 1. The beams are assumed 

to be identical, each having the length L, the cross-sectional 

area A, the Youngs modulas E and the moments of intertia and 

I around the Y and Z axis (principal axes) , respectively, 
z 

The deformation of each joint is described by the displacements 

u, v and w in the X^, ^3 directions, respectivelv and by 

the rotations 9 , , and 9 around the axis X, , X and 

xl x2 x3 123 

respectively. Here the rotations are considered to be positive 
in the counterclockwise direction. 

Using symmetry arguements reveal that this model is orthotropic 
and that a 90° rotation in its plane will not alter its behavior 
[11]. These conditions reduce its general stress-strain relations 
( 1) to 
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which has the four independent constants ^1122' ^1313 


and C 


1212 * 


The actual values of theije constants arc derived 


usinq the direct method of analysis. This consists of isolatino the 
representative repeating cell, figure 2a, loading it at its 
nodes and equating the displacements of these nodes to the 
displacements of the edges of the equivalent continuum plate under 
the same loading conditions. The appropriate loading conditions 
for calculating and shown in figures 2b and 2c, 

those pertaining to calculating C^2i2 shown in figure 2d and 

2e and finally those used in calculating ^^213 shown in 

figures 2 f , 2 g , 2 h . In the first and second loading conditions, 
we are dealing only with the "in-plane" displacements of the 
lattice; while in the third loading condition we are calculating 
the "off-plane" displacement. 

Since each member is shared by two neighboring cells, 
its effective cross sectional area an.ii moments of inertia must 
be half of the corresponding values in the original lattice. 

Under the present loading conditions, matrix structural methods [ 3 ] are 
utilized to solve for the displacements and rotations of each 





individual node. Specifically for fi^gure 2b, we obtain 


Uj . U3 


- V3 . - 


( 8 .b) 


and from figure 2 d, we get 
^ • U2 0 ^ • ^^2 


(9. a) 


u . u » IkL. 
3 4 ^lEI^ 

z 


(9.b) 


Similarly, the displacement in figure 2h is found to be 

^2 “ ^3 “ " PL^/(3EIy/2) (10) 

Figures 2.c, 2.e and 2.f display the equivalent sauare continuum 
element of side length L and thickness h subjected to normal 
stresses, in-plane shearing stresses, off-plane 

shearing stresses, "^^ 3 / respectively. The displacements of the 
plate element due to the normal stress are 


0 -L V 
2 e 


1 E 


and the one due to the in-plane shearing stress 


Ti2 L 


3 G. 


while the displacement due the off-plane shearing stress is 


given as 


c _ _ 
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where E is the effective modulus of elasticity of the equivalent 
e 

orthotropic continuum in the X ^ and the X ^ direction, is the 

efffictive Poisson's ratio of the continuum between the X and 

direction, Gj ^2 in-plane shear modulus and is the off- 

plane shear modulus. The relations between ^ of equation (7) 

and E , V , G,> and G, , are 
e e 12 13 




1122 

'1111 


(14 .a) 


^12 " ^1212 ' ^13 * '^1313 


(14. b) 


By equating the displacements of the plate element with the 
corresponding displacements of the representative unit cell while 
insuring that the total foice on the unit cell equals the total 
force on the plate element for each loading condition yields 


^1111 


AE 

lH 


^1122 = ° 


(15. a) 


^1212 

^1313 


6EI 


3EI 


i7h 


1 


(15. b) 
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b) (0, ^ 60° layups) 

For the (0, ♦ 60°) layup of figure 3, we shall assume that 
all members are identical and have the same geometrical and 
material properties L, A. 1^,1^ and E. The isotropic nature of 
the (0, ^ 60°) configuration (see [2,11]) dictates additional 
restrictiors on the stiffnesses coefficients of the equivalent 
continuum. The appropriate oroperty matri.y is 



r- 





“■ 



1 r 

^1111 

^1122 

0 

0 

0 

0 



^1122 

^1111 

0 

0 

0 

0 



0 

0 

0 

0 

0 

0 



0 

0 

0 

c 

0 

0 

(16) 





1313 





0 

0 

0 

0 

^1313 

0 



0 

0 

0 

0 

0 

1 ^^llll'^1122^ 


which has the three 

independent constants 

^1111' 

Cii 22 and 

The actual values of 

these 

constants 

are 

derived 

using the same 


method outlined above. 

The appropriate loading conditions for calculating and 

^1122 shown in figures 4a and 4b, and those used in calculating 

^1313 shown in figures 4e and 4f. The representative unit 

cell for this layup is shown figure, 4. a. Since the diagonal 
members are shared by two neighboring cells ^ their effective 
cross sectional properties are half those of the chord member. 
v;ith these loading conditions, matrix structural methods are 
utilized ..again to solve for the displacements of each j.ndividual node. 
Specifically, from figure 4c we obtain 
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(17. a) 


(17. b) 


(18) 


Figure 4.b and 4.£ display the equivalent rectangular continuum 
element of side dimensions L x and thickness h, subjected to 

normal streses and off-plane shearina stresses, respectively. 

The displacements of the plate element due to the normal stress 


are given by 


a ^ L 


5 ^ = - /I 


and the displacement due to the off-plane shearing stress r 
is given as 


5 = - 


Using the relations between C. and E , v and G. ^ as given in 
(14) equating the displacements of the plate element with the 
corresponding displacements of the representative unit cell and 
insuring that the 'total force on the unit cell equals the total 
force on the plate element for each loading condition yields 
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(21. a) 


^1111 


3.TEA . 

+ ^ 

4 Lh L-^h 


^1122 


A ea 


(21. b) 


3/3 El 

‘=1313 “ 


(21. c) 


-i2 


IV . Applications 

In this section wo present applications to our con ^.ruction 
procedure as outlined in Sections II and III. T^le models which 
we shall discuss constitute two-dimensional and three dimensional 
beam-like structures, respectively. 

a) Two-Dimensional Structures ; The (0°, 90°, + 45°) layup 

The (0°, 90°, ^ 45°) grid shown in figure 5 is constructed 
from two basic square grids inclined at an angle of 45° and having 
the geometrical properties L, E, A, I , and I and Lk'T, E , , A,, 

y z Cl Q 

^yd' ^zd' respectively. 

The four independent constants for the first (i.e., 0°, 90°) 
basic square grid with respect to its local system of axis are 
given in (15) ; while those corresponding to the ^ 45° square grid 
with respect to its own local system of axis are 


(C 


1111 ^ ^ 

«- 


E, A 
d 

l.T 


d 

h 


^^1122 



(22. a) 


^^ 1212 ^ 


3 E , I , 
d zd 

.1 h 


(22.b) 


(c ) = — 

2 2/2’ h 

The direction cosines of the local system of axis of 
with respect to the fixed coordinate system of axis 
are defined according to (4) as 


(22. c) 

the +45° grid 
'^1' 
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Xi X2 X 3 

(x 1 ) 

r 

r 

0 

^ 2 

/5 /5 

(x 5) 

- i- 0 

^ 2 

/2 /5 

(X 

2 

0 0 1 


Substituting from (15), (22), and (23) into (5) and 
results yield the final properties of the 0°, 90°, 


= EA ^d^d ^ ^^d^zd 
1111 lH 2/5 Lh /2 L^h 


^d^d ^^d^zd 


'1122 


2 vl Lh /2 L h 


^d^d 


6 EI 


1212 


2r'T Lh L h 


L^h 2/5 L^h 
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( 23 ) 


summing the 
45° layup as 

(24. a) 
(24. b) 
(24.0 
(24. d) 



b) 


Three-Dimensional Structures: 


(Octetruss Structures) 


The smallest generating (repeating) unit cell of the octetruss 
structure is shown in figure 6. It is a di amond- like element with 
each of its sides having tho length L and being shared by two 
neighboring cells. The octetruss structure is shown in figure 7 
with respect to the coordinate system arrangement shown in figure 8. 
For further details of the geometric characteristics of this kind 
of structure the reader is referred to [1,2]. In the present 
analysis, the octetruss structure is considered to be composed of 
"beam elements." Examination of this structure reveals ^hat it can 
be construct 2 d from the superposition of different planes. Specific- 
ally, it can be constructed from the three repeating sets of (0°, 90°) 
basic planar grids having different orientation in space, as shown 
in figure 9. The stiffness coefficients for each of the (0°, 90°) 
basic grid with respect to its local system of axis are given in 
(15) where h now stands for the distance between the parallel 
(0°, 90°; Layers; its value is thus given by 


h = <251 

. 2 

The direction cosines of the local system of axis of the three 
basic (0°, 90°) planes with respect to the global system of the axis 
of figure 8 are defined according to equation 4 as ~ 1,2,3 

by 


1/2 

1/2»T* 

* 2/3 

-1/2 

‘ 3/2 

0 

-1/,I 

-»T76 

»T73 


(26a) 
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(6ij)2 
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0 
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✓T 
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-.771 

.T7T. 

-/273 


(26b) 


/T73 


:6c) 


Substituting from (26) into (5), using (15) and summing according 
to (6) yields the final properties of tne octetruss structure 
with respect to coordinates of figure 8 as 

m 

C 
C 
C 
C 


(C 


ijkl 


1111 
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c 
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Sl23 

0 

0 

1122 

^1111 
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^2323 

0 

0 

0 

0 

0 

0 
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where 
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(28. a) 


(28. b) 


(28. c) 


(23. d) 


( 28. e) 
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■1212 ' n~ 7? ^ "T^ ^ T 
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(28. f) 
(28. g) 


Notice that (27) constitutes a modification of our previously 
reported result in (1) which are reflected in the appearance of 
the bending rigidties of the members. Notice also that there is 
no change in the number of the independent constants which can 
also be deduced from symmetry (1,2). Examination of the results 
(28) indicates that “ ^^1111 ” ^1122^'^^ hence the octetruss 

is transversely isotropic, as is expected. 

Remar)< 

By reexamining figure 7 we can see that the same structure can 

also be constructed from four different repeating sets of (0°, + 60°) 

basic planar grids. In this case of construction, each member will 

be shared by two different basic grids. Since I and I are the 

y z 

moments of inertia of the cross section of the beam around two 
principal axes and since each beam is shared by two different basic 
grids, we must have two sets of principal axis for each cross 
section; this can only sense for circular cross-sections. Thus, 
constructing the properties of the octetruss from those pertaining 
to four (0°, + 60°) layup 15 restrictive in that only beams with 
circular cross-section can be treated. This was actually done and 
its results were found identical to (27) and (28) when the later 
are also soecialized to I = I . 
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b - Equivalent Plate Sign Convention. 

Sign Convention for the Displacements in 
the Equivalent Continuum Plate Model. 


Sketch 1^ 



Figxire 1. The (0®,90*) Lattice 
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Figure 2. The Loading Conditions of the Representative Repeating 
Cell for the (0®,90*) Layup used to detenaine the 
Stiffnesses Coefficients of the Equivalent Continuum, 






Figure 3. ( 0*, + 60* ) Layup . 
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Figure 4. 


The Loading Conditions of the Representative 
Repeating Cell for the (0“,^60*) layup used 
to determine the Stiffnesses Coefficients of 
the Equivalent Continuum. 



Flgur« 4 (cont.) The Loadim^ Conditions of the Representative 

I 

Repeating Cell for the (0*,^60*) Layup used to 
determine the Stiffnesses Coefficients of the 
Equivalent Continuum . 




rigura 5. Tha (0»,45*, * 90*) Lattica . 



Figur« 6. Smallsst Repaating Element of the 
Octetruss Struct\ire . 
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Figure 7. Three-Dimensional Octetruss Structure rlewed with 
respect to the Coordinate System of Figure 8. 



Figure 8 . Direction Cosines of the Octetruss . 
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First Grid 
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Third Grid 
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